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SENSITIVITY ANALYSIS OF NEURONAL DYNAMICS BASED ON 

ADDITIVE STDP RULE 

SUBHAJIT SENGUPTA*, KARTHIK S. GURUMOORTHY*, AND ARUNAVA BANERJEE * 

Abstract. Spike Timing Dependent Plasticity (STDP) is a Hebbian like synaptic learning rule. 
The basis of STDP has strong experimental evidences and it depends on precise input and output 
spike timings. In this paper we show that under biologically plausible spiking regime, slight variability 
in the spike timing leads to drastically different evolution of synaptic weights when its dynamics are 
governed by the additive STDP rule. 
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1. Introduction. Spike Timing Dependent Plasticity (STDP) [15, 2, 17, 11, 4, 6, 
5] is a temporally asymmetric form of Hebbian learning which has been experimentally 
observed in a wide spectrum of species. It is widely believed that STDP plays a 
key role in learning [1, 7], information storage [10], coincidence detection [9] as well 
as the neuronal circuits development in the brain. STDP rule also has attractive 
features such as synaptic competition and rate stabilization [15]. It is induced by the 
strong temporal correlation between the spike timings of presynaptic and postsynaptic 
neurons. The synaptic modihcation depends on the precise timing and the order of 
the input and output spikes. The synaptic connection between two neurons is more 
likely to strengthen if the presynaptic neuron fires shortly before the postsynaptic 
neuron and similarly weakens if the inverse situation happens [13]. As STDP depends 
on specihc spike timing, it would be interesting to investigate a scenario where those 
spike timings are perturbed slightly. In other words, we seek to answer the question: 
“Is STDP robust enough to make synaptic weights insensitive to slight variations in 
input (or output) spike timings?” Our analysis reveals the contrary where small 
initial perturbations have drastic effects on synaptic weights when its dynamics are 
governed by additive STDP rule. We would like to accentuate that though the fragility 
of the additive STDP rule concerned with the stabilization issues has been previously 
reported in the literature [16, 17, 14], our analysis is different and directed from a 
dynamical system viewpoint. 

In order to investigate this, we studied the dynamics of a single output neuron with 
multiple synapses whose weight changes are dictated by the STDP rule. At the initial 
STDP time-window we give a tiny perturbation on input or output or both spikes. 
Later we do not allow any other perturbations in this system and let the system evolve 
in time. Our formal study shows that under biologically plausible spiking regimes the 
evolution of synaptic weights would be very different in this perturbed case than the 
original unperturbed system due to this initial small perturbations. For our analysis, 
we framed an equivalent problem to this by fixing the spike timings in the initial STDP 
time-window and adding a small perturbation on the synaptic weight itself. We argue 
that this two situations are exactly same because due to those small perturbations in 
input or output spike timings, the STDP update for synaptic weight would be slightly 
different in both cases, which will serve as the initial synaptic weight perturbation. 
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In the entire analysis, we work with the perturbation vector associated with synaptic 
weight and prove our results. Our analysis considers a spiking neuron as an abstract 
dynamical system as described in [3]. In [3] the dynamical system was studied in the 
context of recurrent spiking neuron network assuming a fixed synaptic weight. In our 
case we study the dynamics of a single neuron particularly focusing on how STDP rule 
on excitatory synapses affects initial perturbations associated with synaptic weights. 
We study the effects of the six events, namely (1) birth of an output spike, (2) birth 
of an input spike, (3) death of an output spike, (4) death of an input spike, (5) eeiling 
and (6) floor. The last two events are the manifestation of the additive STDP rule. 


1.1. A simple experiment. We start with the following simple experiment. 
Consider a single output neuron innervated by 1000 input synapses among which 
20% are inhibitory. Spiking dynamics in the system are governed by the first-order 
Spike Response Model [8]. At each timestep, the output neurons membrane potential 
is computed as the sum of all Post-Synaptic Potentials (PSP) of its input spikes, 
summed with the After Hyper-Polarizing (AHP) effects of previous output spikes. 

_ 

Form of PSP is given by Pflxl^) = —where xl^ denotes the elapsed 

dy 

time at the synapse since the occurrence of the event kj (the notion of an event is 
made clear below), d(> 0) denotes the distance (in dimensionless units) of the synapse 
from the soma, and /3(> 0) and 7 p(> 0) control the rate of rise and fall of the PSP. 
The after hyper-polarizing potential (AHP) effects are determined completely by the 


time since the departure of each output spike. The AHP is given by ) = Re 

where R{< 0) denotes the instantaneous fall in potential after a spike and 7a(> 0) 
controls its rate of recovery. represents the elapsed time since the event kj took 
place. In all experiments, we set /3 = 6 (dimensionless), r = 15msec, R — — lOOOmH 
and 7 a = 1.6msec. For additive STDP, A+ and A_ determine the maximum synaptic 
update amounts, and and t_ determine the steepness of the exponential response 
curves. We use Ap = 0.01, A_ = 0.0105 and t_|_ = t_ = 20msec. The synaptic 
weights between inputs and output are initialized randomly between Wmax and Wmin 
whose values are 40 and 0 respectively. 


We simulate two runs of the same system with identical Poisson spike train in¬ 
puts. We run them simultaneously till the synaptic weights reach to a steady state 
distribution. Once the weights have stabilized, we slightly perturb the weight of a 
single synapse in one of the system. In the Figure 1.1a, we show the 2D plot of the 
synaptic weights for the perturbed and the unperturbed system in the x-axis and 
y-axis respectively and they are same except for a single synapse (the red dot shifted 
away from the diagonal corresponds to the perturbed synapse). We then let both the 
perturbed and the unperturbed systems to evolve with the synaptic weights governed 
by the additive STDP rule. After some long enough time, final synaptic weights for 
both the systems are noted and plotted in the Figure 1.1b. Equal synaptic weight 
in both the systems are shown in ‘blue’ and synapses with different weights in both 
the systems are shown in ‘red’. We detect that it has been deviated from the x = y 
line for almost all the places - which clearly indicates sensitivity to initial condition. 
The claim is further strengthened by the plot in the Figure 1.1c, where we show 
the normalized sum of squared difference of synaptic weights between perturbed and 
unperturbed systems and find that to constantly increase. 
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(c) Normalized sum of 

(a) Both systems - Initial (b) Both systems - Final squared differences 


Figure 1.1 


1.2. Brief exposition of our result. We now provide a brief exposition of 
why the dynamics governed by the STDP update rule might be sensitive to small 
perturbations to synaptic weights. The formal analysis is described in detail in the 
later sections. Let us consider the two events - hirth of an output spike followed 
by a birth of an input spike to happen after a tiny increment of weight in one of 
the input synapses. Due to larger weight in the perturbed synapse, the perturbed 
neuron will spike little earlier compared to the unperturbed one. Hence the time 
difference between the current output spike and the past input spikes is smaller in 
the perturbed system giving rise to larger increment in the synaptic weights via the 
additive STDP rule. This enhances the differences in the synaptic weights. Now for 
an impending input spike, the most recent output spike in the perturbed system is 
moved little earlier to the past (as it spiked earlier). Hence the negative STDP update 
would be less in the perturbed system in comparison to the unperturbed one, thereby 
further enlarging the weight perturbation. This gives us an inkling that the additive 
STDP rule might be sensitive to small weight perturbations, viz. sensitive to small 
perturbations of input or output spike timings. 

2. Phase space dynamics. Before we analyze the aforementioned events in 
greater detail, let us consider phase space of our dynamical system. In order to 
track the weight value at the synapses, we need to consider all the past input and 
output spikes within the neuron’s synaptic efficacy window. As an input spike picks 
up the current weight associated with the input synapse at which it took birth, it 
is convenient to associate past synaptic weights with its imminent input spike. The 
phase-space of our dynamical system will thus consist of all past weights associated 
with the input spikes, current synaptic weights and the past output spikes within the 
synaptic efficacy window. 

As the number of input and output spikes in the synaptic efficacy window gen¬ 
erally varies with time, so will be the dimensionality of our phase space. Since the 
maximum number of spikes in a synapse that any neuron can have within its efficacy 
window T is bounded by L = [y], where r is the refractory period for that neuron, 
an easier fix to the changing dimensionality problem is to have our phase-space be of 
constant dimension D = m + mL + L, where m is the number of input synapses. The 
quantity mL and L are for the input and output spikes respectively. If a synapse has 
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less than L spikes, the remaining entries are set to 0. 

We study the effects of the six events, namely (1) birth of an output spike, (2) 
birth of an input spike, (3) death of an output spike, (4) death of an input spike, (5) 
ceiling and (6) floor. The last two events are the manifestation of the additive STDP 
rule. We will discuss effect of each of these events in the following subsections. 

2.1. Birth of an Output Spike. Let T be the threshold for the membrane 
potential and let kj correspond to the event when the spike occurred at either 
the input synapse or the output synapse. As before, let and Xq^ denote the 
time of occurrence of the spike in the input synapse and the output synapse 
respectively and let denote the weight associated with the input spike x^^ with 
values between Wmax and Wmin- Let Ui and no denote the number of spikes in the 
input synapse and output synapse respectively within the spike efficacy window. 
Consider the event when the unperturbed system produced an output spike where we 
have 

m rii no 

(2.1) ^ ^ ^ ) + E Po(4n- 

i=l j=l j = l 

Let Awf denote the initial weight perturbation for the synapse just before this 
event. Let Aa;§ denote the perturbation in the output spike time due to this weight 

A) * 

perturbation. Let the perturbations associated with the past output spikes Xq be 

k ■ 

denoted by Axq^ . We then have, 

m rii no 

(2.2) T = E E^'^^ + + E Po{xq^ + Axq^ — Axq). 

i=i j=i 

Using I'’* order Taylor expansion, equating with (2.1) and ignoring any second order 
term we have. 




dP, fc, , ^ aPo 


i=l j=l 


dr 

= 1 '^■^0 


-W, 


E 


r\ kj 'b r\ kj 

i=ij=idxfl j=i oV 


(2.3) 

Hence the resulting output perturbation because of changes in the synaptic weight is 
given by, 


(2.4) 

where. 


A4 = 


EIli YTjU P.(^r) + YTU 

OXf) 


(2.5) 


a = EE 


fc, , aPo 


-W, 


E 


r\ kj 'b r\ kj 

i=ij=idxfl i=iOV 


From the additive STDP update, an output spike increases the synapse’s weight 
by 

= Wmax where A_|_, t+ are constants for positive update of additive 

STDP. The new updated weight for synapse after this event equals . 
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Now consider the perturbed system. Let Wi^ denote the perturbed weight at the 
synapse before the occurrence of this output spike, i.e, Wi^ = ref + Arcf. Let 
denote the new weight in the perturbed system after the event. Then from 
STDP rule we have, 

(2.6) wT+i = Wj^ + Wniax A+e "+ 

i=i 

Using the 1®* order Taylor expansion and plugging the value of Wi^ we get, 

rii _J^3 . j. 

(2.7) Arcf+1 = ^7- w’l+^ = Au;f + Wma. Y^A+e^^^. 


Substituting the value of Axg from Equation 2.4 we can write 


Au;f+i = Aa;f + ( 




T+ 


(TZi e;=i a^^p,(4o + yTU 

A 


( 2 . 8 ) 


In order to analyze the sensitivity of the STDP rule, it suffices to study the dy¬ 
namics of the perturbation vector during the occurrence of these events. As briefly 
mentioned before, the perturbation vector has three components - (i) weight pertur- 

^ (fc) 

bations associated with m synapses (Awg ) (ii) weight perturbations associated with 

— (k) -* (k) 

all the input spikes {Awj ) (iii) output spikes perturbations(Aa;Q ). Figure 2.1 shows 
all the components with appropriate dimensions. 



Figure 2.1: Components of the perturbation vector 


Let 


(2.9) 




T+ 


i=i 


and = and 


1 dPa 


Zn = -r 


^ dxn 


Define the vectors 


u^^^ = [ulul- 


= [zo ^^, zo ’^^, • • •, zo*^-o , 0,0, • • • 0], 


= It! 


1 5 


^ni ^ fci ^^2 r\ 

■^Vl , 0 , •••,^^2 


, Vm , 0, • • •] 


and the perturbation vectors 

= [Awf,- • •, Aw^], Axo^^^ = [Aa:o''\ • • •, Aa;o'‘"o,0,0, • • • ,0] and 
Awf'’ = [Au>i 1, • • •, AwY , 0, • • •, Aw 2 ^ , • • •, AwY , 0, • • •, Aw ^, • • •, Awm ’", 0, • • •]. 
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Note that in the above definition, some of the vectors have zeros appended to it in 
order to maintain the constant dimensionality of the phase-space as described under 
Section 2. Weight update takes place only at the synapses while the weights associated 
with input synapses remain intact. The new output perturbation Axq gets added to 

(fe) 

the output perturbation vector Axq . The matrix governing the dynamics of the 
perturbation vector is shown in Figure 2.2 where / represents the identity matrix of 
appropriate dimension and Ji is an identity matrix appended with a last column of 
all zeros. 

Firstly note that is positive as it is sum of exponential each multiplied by a 
positive constant A+. For excitatory synapses, the numerator term in the definition 
of namely, Pi{xij), is positive as it is just the post synaptic potential (PSP) 
corresponding to the input spike Xij. For similar reasons, the numerator will be 
negative if the synapse is an inhibitory one. Since STDP rule doesn’t apply for 
inhibitory synapses (as their synaptic strength doesn’t change) we may set = 0,Vj 
as its corresponding Aw^ ’' = 0 - the latter is true as there is no perturbation at an 
inhibitory synapse before the event. In other words, as Aw^^ = 0 independent 
of uf”, without loss of generality we can set =0. Additionally uf = 0 for an 
inhibitory synapse. So we can restrict ourselves only to excitatory synapses in our 
subsequent analyses. 

k * 

To see > 0 for an excitatory synapse, observe that its denominator A is the 
derivative of the membrane potential at the time of the output spike in the unper¬ 
turbed system and hence is positive. Also the components of the namely Zq^ are 
strictly positive as its numerator is just the derivative of the AHP which is a strictly 
increasing function. Hence we notice that all the entries in the perturbation matrix 
are non-negative, a fact that is used in Section 3 to show the sensitivity of the additive 
STDP rule. 

2.2. Birth of an Input Spike. For our analysis we consider the event birth of 
a input spike to have occurred after a set of past output spikes in the neuron’s STDP 
efficacy window. Hence an input spike will decrease the synaptic weight at the p*^ 
synapse where it took birth. The decrement based on additive STDP rule is given 

n ^0 

by Cp = '^maa: where denotes the past output spike timings, A- 

and r_ are positive constants for negative STDP update. Closely following the steps 
delineated in Section 2.1 and using the 1®* order Taylor expansion, the synaptic weight 
difference at the synapse between the perturbed and the unperturbed systems after 
the arrival of the current input spike can be computed as 

(2.10) Awp^^^ = Awp'" Wraax 6 ^ ° ■ 

i = l 

Moreover, as the newly born input spike will pick up the previous synaptic weight 
and its corresponding weight perturbation from the synapse, a new entry needs to 
be added to Awi at the (p — 1)L -I-1 location with a value of Awp. The perturbation 
matrix is given in Figure 2.3. J 2 is a zero matrix except having 1 at the ((p— l)T-f l,p) 
entry for addition of Awp element to Awi at the (p — 1)L-|-1 location. Both /i and I 2 
are identity matrices of dimension (p—1)L and {m—p)L respectively corresponding to 
unaltered weight perturbations associated with input spikes in other synapses. K is 
another identity matrix appended with a first row and last column of all zeros. As all 
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added 

at ((p — 1) * Z, + 1)'^ entry 



deleted entiy 


Figure 2.2: Matrix for birth of an output Figure 2.3: Matrix for birth of an input 
spike spike 


the synaptic weights except for the synapse remain unscathed, all rows in C barring 

the row are zero and the entries in the p*^ row are given by Cp. = A_ e . 
Importantly notice that all entries in the perturbation matrix are again non-negative. 
We will make use of this fact in our subsequent analysis. 


2.3. Death of an output and input spike. These events do not change the 

(fe) 

synaptic weights. Only the entries in the Uq^ and {L{p— 1) +npY^ locations in Aa:o 
(fe) 

and Awj are set to 0 respectively, p denotes the synapse where the event death of 
an input spike took place. The perturbation matrices for these events are shown in 
Figures 2.4 and 2.5. Here, I denotes the identity matrix of appropriate dimension. J 3 
and J 4 are very similar to the identity matrix except they have one row of all zeros 

(fe) (fe) 

for setting the corresponding entries to 0 in Axq and Awj . 



Figure 2.4: Matrix for death of an output 
spike 


Figure 2.5: Matrix for death of an input 
spike 


2.4. Ceiling and Floor Events. It is well known in literature [15] that when 
the dynamics of the synaptic weights are governed by the additive STDP rule, an 
artificial cap needs to be imposed on the maximum (minimum) values that the weights 
can take. During the occurrence of these ceiling (floor) events corresponding to the 
weights attaining the maximum (minimum) values, the corresponding STDP rule 
is not applied and the weights values stays put. As a result, the entries of Aws 




















































































































































































































































corresponding to the synapses which faces these events will be set to zero in the worst 
case nullifying the perturbation. 


3. Sensitivity Analysis. We begin our sensitivity analysis by demonstrating 
that except for the ceiling and the floor events, the additive STDP rule is sensitive to 
small synaptic weight perturbations. Since the occurrence of the ceiling and the floor 
events are stochastic in nature and completely depends on the input spike distribution, 
the neutralizing effect of these events on the perturbation vector needs to be modeled 
via a stochastic process as described in the subsequent section. The overall sensitivity 
of the STDP rule will depend on the interplay between the enhancing effects of the 
birth and death events and the squashing effect of the ceiling and floor events on the 
perturbation vector. 


Let denote the perturbation matrix for the event and let represent 

the perturbation vector after the event, i.e., = It is worth 

recalling that all entries of A^^'> are non-negative irrespective of the type of event. 
Additionally the perturbation vector is made of three components namely = 

r A ^ (^) A (^) 

Awj Ax \ as shown in Figure 2.1. 

Let = A^^'> * • • • * denote the combined effect of these perturbation 

matrices after the occurrence of n events. In order to prove the sensitivity of the 
additive STDP to small perturbation, it suffices to show that as n —>■ oo, the Frobe- 
nius norm of the matrix goes to infinity. In order to show this, without loss 
of generality we take a specific initial perturbation vector y*-^) given by y^^^ = 
{O, • • •, 0, Awp, 0, • • •, 0} , 0, 0 where only the synapse (for any p) has a non-zero 


positive perturbation Awp. Both the input weight perturbation and the 

output perturbation vectors are set to 0. Then ||y ^^^||2 = HAwg^ ^||2 = Awp. 

We now show that this choice of y^^^ works, i.e., Hy^"'*'^^II 2 > HAws^ ^||2 > 
llAics^ ^||2 = ||y*'^^|| 2 , for sufficiently large n. Specifically we prove that the norm 

(k) 

of Awg always increases with the occurrence of specific type of events and is never 
decreased by any of these events (except for the ceiling and floor events for which we 
gave the stochastic analysis in the next section). 


Since all the entries of A^^'> are non-negative^ it is easy to see that all entries of 
^(fc+i) are also non-negative even after the occurrence of any of these events. Focusing 

on the dynamics of the first component of namely Aws , firstly, observe 

that the occurrence of death of either an output spike or an input spike doesn’t affect 
. - (fc+i) 

Aws as the synaptic weights remain unruffled. Secondly the birth of the output 

or input spike always increases the norm of Aws as either all its components (in 
case of birth of an output spike) or one of its components (in case of birth of an input 
spike) gets a positive increment. This is easily verifiable from Equations 2.8 and 2.10 
because of positivity of y^^^. In fact we can claim that this positive increments are 

-» (fe) 

at least e which are bounded away from 0. In short, the ||Art;s II 2 increases with 
the occurrence of birth of either the output or input spikes and doesn’t change with 
the death of these spikes. Hence for sufficiently large n with 

increments bounded away from 0, which in turn suggests that the Frobenius norm of 
Afo) will increase and converge towards infinity. 
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3.1. Characterizing the Underlying Stochastic Process. The aforemen¬ 
tioned analysis requires an investigator who has access to the system of this spiking 
neuron from which he draws data to compute the perturbation matrix for all 
the birth, death, ceiling and floor events. With the absence of the last two events, 
our neural system is sensitive to small weight perturbations irrespective of what¬ 
ever statistics it might have. The occurrence of the ceiling and floor events and its 
stochastic nature warrants the development of a stochastic model for our analysis. We 
assume that both the input and output spike distribution is a homogeneous Poisson 
process with the mean rate of Ai„ and Aon* respectively. Let pi, p 2 , Ps, P 4 be the 
steady state probability of occurrence of birth of an output spike, birth of an input 
spike, death of an output spike and death of an input spike respectively. Since the 
occurrence of the ceiling event requires that the synaptic weights attain or exceeds its 
maximum value, it can take birth only after the birth of an output spike. Similarly, 
the floor event can occur only after the birth of an input spike. Let 71 and 72 be the 
probability of occurrence of ceiling and floor events given the occurrence of the birth 
of an output and input spike respectively and also assume that with probability 73 , 
a synapse individually experiences the ceiling event given that the ceiling event has 
already occurred at that time point. We now describe the stochastic process for the 
generation of the perturbation matrices corresponding to each of these events. 
The process is begun by choosing integers ni{i = 1, ■ ■ ■ ,m) and no randomly over the 
range [0, L] to denote number of input spikes in the synapse and output spikes 
inside the synaptic efficacy window respectively. 

In case of birth of an output spike, the perturbation matrix is shown in 
Figure 2.2. The entries of the matrix are functions of u*, and Zq^ given in Equa¬ 
tion 2.9 which in turn depends on the input and output spikes-through the PSP and 
AHP functions and its derivatives-and the weights associated with the input spike 
train. Hence the matrix can be populated by sampling from the input and output 
spike times which follow Poisson distribution and from the synaptic weight distri¬ 
bution which is known to be bimodal [15] and then computing the and Zq^ 

values. Care should be taken to ensure that the value of A given in Equation 2.5 is 
positive. We also assume that the choice of birth and death events are independent 
of the entries of the matrix. In other words the distribution of the entries obtained as 
function of input and output spikes and synaptic weights is independent of the birth 
and death events. If the ceiling event follows this output event (with the probability 
of 7 i), then all the entries of certain rows in the constructed perturbation matrices 
will be set to zero. 

In case of birth of an input spike, the choice of the synaptic index p at which 
the event is expected to occur is randomly chosen between 1 to m (m is the number 
of input synapses) independent of the previous events. Once the choice is made, the 
corresponding perturbation matrix A^^'^ shown in Figure 2.3 is populated by sampling 
from the output spike time Poisson distribution and then computing the Cp. values. 
The distribution of Cp. values obtained as a function of output spike train distribution 
is assumed to be independent of the past events. If the floor event follows this input 
event (with probability 72 ), then all entries of row of A^^'> are set to zero. 

In case of a death of an output event, the perturbation matrix is an identity 
matrix except that the entries of the row corresponding to the last existing output 
spike is set to zero. Similarly, in case of a death of an input event, the synaptic index 
p at which the death event is expected to occur is again chosen randomly between 1 
to TO. The perturbation matrix is an identity matrix except that the entries of the 
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TImesteps 


Figure 3.1: Frobenius norm of A^"^. (i) Blue: High ceiling and floor events, (ii) Red: 
Low ceiling and floor events 


row corresponding to the last existing spike in the synapse is set to zero. 

In Figure 3.1 on left, two very different scenarios of simulation are presented. 
Both curves indicate the Frobenius norm of computed in various timesteps. The 
blue curve shows the rise of the Frobenius norm, which clearly reaffirms our analysis 
as we have lesser probabilities for ceiling and floor events (71 = 72 = 0 . 1 ) and 
lesser probability for ceiling event for each synapse (73 = 0.1) as well. On the other 
hand, the red curve shows minimum variance with a decreasing trend as we have 
7 i = 72 = 73 = 0.9, which decreases the Frobenius norm by generating more and 
more ceiling and floor events. 

4. Discussion. We studied the dynamics of the synaptic weights of a single neu¬ 
ron where the weights are modified according to the additive STDP rule. Interestingly, 
we noticed that the individual synaptic weight values between the perturbed and the 
unperturbed system tend to vary significantly as they evolve in time demonstrating 
the sensitivity of the additive STDP rule to tiny perturbations. If either the synaptic 
weight, or the input and output spike times of a neuron were to be infinitesimally 
perturbed at a given point in time, additive STDP would cause the neuron to reach 
vastly different states, even when the subsequent input spike train is identical. We 
also noticed that as the frequency of the ceiling and floor events are increased (when 
the weights reach their maximum or minimum values respectively)—scenarios where 
the STDP rule is not applied—the system is more or less stable and becomes insen¬ 
sitive to perturbations which is quite counter-intuitive. In the future, we would like 
to investigate other kinds of STDP rule like multiplicative, add-multiplicative, mixed 
STDP rules. We would like to recover the criterion which decides the sensitivity of 
the system. In a seminal article [12], it has been shown how one can convert a chaotic 
attractor to any one of a large number of possible attracting time-periodic motions 
by making only small time-dependent perturbations of an available system param¬ 
eter. Another fruitful research prospect would be to study how to counter-balance 
sensitivity effects and drive the system to the desired state. 
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